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00 ' The eikonal description of breakup reactions diverges because of the Coulomb 

■ interaction between the projectile and the target. This divergence is due to the 
. adiabatic, or sudden, approximation usually made, which is incompatible with the 

infinite range of the Coulomb interaction. A correction for this divergence is analysed 
by comparison with the Dynamical Eikonal Approximation, which is derived with- 
out the adiabatic approximation. The correction consists in replacing the first-order 
term of the eikonal Coulomb phase by the first-order of the perturbation theory. This 
allows taking into account both nuclear and Coulomb interactions on the same foot- 
ing within the computationally efficient eikonal model. Excellent results are found 
for the dissociation of ^^Be on lead at 69 MeV/nucleon. This Coulomb Corrected 
O ' Eikonal approximation provides a competitive alternative to more elaborate reaction 

^ ■ models for investigating breakup of three-body projectiles at intermediate and high 

energies. 

■ PACS numbers: 24.10.-i, 25.60.Gc, 03.65.Nk, 27.20. +n 
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I. INTRODUCTION 



Halo nuclei are among the most peculiar quantum structures [H, 0, [sj . These light neutron- 
rich nuclei exhibit a very large matter radius when compared to their isobars. This extended 
matter distribution is due to the weak binding of one or two valence neutrons. Thanks to 
I their low separation energy, these neutrons tunnel far inside the classically forbidden region, 
and have a high probability of presence at a large distance from the other nucleons. In a 
simple point of view, they can be seen as very clusterized systems: a core that contains 
most of the nucleons, and that resembles a usual nucleus, to which one or two neutrons are 
loosely bound, and form a sort of halo around the core j4|. The ^^Be, ^^C, and ^^C isotopes 
are examples of one-neutron halo nuclei. Examples of two-neutron halo nuclei are ^He, ^^Li, 
and ^^Be. In addition to their two-neutron halo, these nuclei also exhibit the Borromean 
property: the three-body system is bound although none of the two-body subsystems is jsf. 

Since their discovery in the mid 80s m, these nuclei have thus been the focus of many 
experimental [H, S| and theoretical 0| studies. Due to their short lifetime, halo nu- 
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clei cannot be studied with usual spectroscopic techniques, and one must resort to indirect 
methods to infer information about their structure. Breakup reactions are among the most 
used methods to study halo nuclei [l^, 11, T2|- In such reactions, the halo dissociates from 
the core through interaction with a target. In order to extract valuable information from 
experimental data one needs an accurate reaction model coupled to a realistic description of 
the projectile. Various techniques have been developed with this aim: perturbation expan- 
sion [H, l3|, adiabatic approximation jl^, eikonal model [l^, 13, coupled channel with 
a discretized continuum (CDCC) [19i, |20|, i21j], numerical resolution of a three-dimensional 
time-dependent Schrodinger equation fTDS El [221. 



H, M, S H i3, and more recently 
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dynamical eikonal approximation (DEA) [28|, |29|, 

Some of these techniques (perturbation expansion, adiabatic approximation, and eikonal 
model) are based on approximations that lead to easy-to-handle models. Their main ad- 
vantage is their relative simplicity in use and interpretation. However, the approximations 
on which they are built usually restrain their validity domain. For example, perturbative 
and adiabatic models are restricted to the sole Coulomb interaction between the projectile 
and the target. The eikonal method on the contrary diverges for that interaction and can 
be used only for reactions on light targets. The adiabatic, or sudden, approximation made 
in the usual eikonal model is responsible for that divergence. It indeed assumes a very brief 
collision time, that is incompatible with the infinite range of the Coulomb interaction. 

The more elaborate models (CDCC, TDSE, and DEA) are not restricted in the choice 
of the projectile-target interaction. However, they lead to complex and time-consuming 
implementations. First calculations were therefore limited to simple descriptions of the 
projectile (i.e. two-body projectiles with local core-halo interactions). Recently, several 
attempts have been made to improve the description of the projectile. For example Summers, 
Nunes, and Thompson have developed an extended version of the CDCC technique, baptized 
XCDCC, in which the description of the halo nucleus includes excitation of the core (3l| . 
Other groups are developing four-body CDCC codes, i.e. a description of the breakup of 
three-body projectiles, with the aim of modeling the dissociation of Borromean nuclei [s^. 



331 ]. These techniques albeit promising, require large computational facilities, and are very 
time-consuming. 

Alternatively one could try to extend the range of simpler descriptions of breakup re- 
actions. Among these descriptions, the eikonal model is of particular interest. It indeed 
allows taking into account, at all orders and on the same footing, both nuclear and Coulomb 
interactions between the projectile and the target. Moreover it gives excellent results for 



nuclear-dominated dissociations [17|, |29|. Its only flaw is the erroneous treatment of the 
Coulomb interaction. A correction to that treatment has been proposed by Margueron, 
Bonaccorso, and Brink 34 1 and developed by Abu- Ibrahim and Suzuki [35|. The basic idea 
of this Coulomb corrected eikonal model (CCE) is to replace the diverging Coulomb eikonal 
phase at first-order by the corresponding first-order of the perturbation theory [36|. The 
latter, being obtained without adiabatic approximation, does not diverge. The CCE is much 
more economical than more elaborate techniques (a gain of a factor 100 in computational 
time can be achieved between this CCE and the DEA). It could therefore constitute a com- 
petitive alternative for simulating the breakup of Borromean nuclei at intermediate and high 
energies. However efficient it seems, this correction has never been compared to any other 
reaction model. 

In this work, we aim at evaluating the validity and analyzing the strengths and weaknesses 
of this correction by comparing it with the DEA. The chosen test cases are the breakup of 
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^-^Be on Pb and C so as to see the significance of tlie correction for both heavy and hght 
targets. The considered energy is around 70 MeV/nucleon. This corresponds to RIKEN 
experiments 11, l3], with which the DEA is in excellent agreement 28, 29 . 

Our paper is organized as follows. In Sec. [TTl we recall the basics of the eikonal description 



of reactions, and detail the Coulomb correction proposed in Refs. |34J . l35l | . The numerical 
aspects of our calculations are summarized in Sec. IIIII The results for ^^Be on Pb are 
detailed in Sec. [TVl while those corresponding to a carbon target are given in Sec. |Vl The 
final section contains our conclusions about this model. 



II. THEORETICAL FRAMEWORK 
A. Eikonal description of breakup reactions 

To describe the breakup of a halo nucleus, we consider the following three-body model. 
The projectile P is made up of a fragment / of mass ruf and charge Zje, initially bound to 
a core c of mass rric and charge ZcC. This two-body projectile is impinging on a target T of 
mass rriT and charge Ztc. The fragment has spin /, while both core and target are assumed 
to be of spin zero. These three bodies are seen as structureless particles. 

The structure of the projectile is described by the internal Hamiltonian 

Ho = ^ + Kf{r), (1) 

where r is the relative coordinate of the fragment to the core, p is the corresponding momen- 
tum, fief = mcmf/mp is the reduced mass of the core-fragment pair (with mp = rric + mf), 
and is the potential describing the core- fragment interaction. This potential includes a 
central part, and a spin-orbit coupling term (see Sec. IIIip . 
In partial wave Ij, the eigenstates of Hq are defined by 

Hoct)ijm{E,r) = E<Pi,^{E,r), (2) 

where E is the energy of the c-f relative motion, and j is the total angular momentum 
resulting from the coupling of the orbital momentum / with the fragment spin /. The 
negative-energy solutions of Eq. ([2]) correspond to the bound states of the projectile. They 
are normed to unity. The positive-energy states describe the broken-up projectile. Their 
radial part Ukij are normalized according to 

Ukij{r) — ^ cos6ijFi{kr) + sm6ijGi{kr), (3) 

r^oo 

where k = \j2^cfE /h^ is the wave number, 5ij is the phase shift at energy E, and Fi and 



Gi are respectively the regular and irregular Coulomb functions [37 , 

The interactions between the projectile constituents and the target are simulated by op- 
tical potentials chosen in the literature (see Sec. IIIII) . Within this framework the description 
of the reaction reduces to the resolution of a three-body Schrodinger equation that reads, in 
the Jacobi set of coordinates illustrated in Fig. [H 
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^ + i/o + Vpt{R, r) 



^{R,r)=ET^{R,r), (4) 
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where R is the coordinate of the projectile center of mass relative to the target, P is the 
corresponding momentum, fi = mpuiT / ijnp + mx) is the projectile-target reduced mass, 
and Et is the total energy. The projectile-target interaction 



VpT{R,r) = V,T 



R-^r]+VfT 
mp ' 



_ rric 
R+—r 

mp 



(5) 



is the sum of the optical potentials Vct and V/t (including Coulomb) that simulate the 
core-target and fragment-target interactions, respectively. The projectile impinging on the 
target is initially bound in the state (^iojomo '^^ energy Eq. We are therefore interested in 
solutions of Eq. (jl]) that behave asymptotically as 



^(R,r) 



(6) 



where Z is the component of R in the incident-beam direction and rj = ZrpZpe^ / [ATreohv) 
is the P-T Sommerfeld parameter (with Zp = Zc + Zf). 

In the eikonal description of reactions, the three-body wave function \& is factorized as 



the product of a plane wave by a new function \1/ [16|, ll7|, [18 



^(il,r) = e*^^^(i?,r). 



(7) 



where K is the wavenumber of the projectile-target relative motion related to the total 
energy Et by 



Ej 



2/z 



+ ^0- 



(8) 



With factorization ([7]), the Schrodinger equation (jl]) reads 



p2 



H0-E0 + Vpt{R, r) 



^{R,r) = 0, 



(9) 
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where v = hK/ ^ is the initial projectile-target relative velocity. The first step in the eikonal 
approximation is to assume the second-order derivative negligible with respect to the 

first-order derivative vPz- The function ^ is indeed expected to vary weakly in R when the 
collision occurs at sufficiently high energy 

0, H This leads to the DEA Schrodinger 

equation [2^, 2^ 



d - 

ihv — ^{b, Z, r) 



[{Ho-Eo) + VpTiR,r)]<il{b, Z,r) 



(10) 



where the dependence of the wave function on the longitudinal Z and transverse b parts 
of the projectile-target coordinate R has been made explicit (see Fig. [1]). This equation 
is mathematically equivalent to a time-dependent Schrodinger equation with straight-line 
trajectories, and can be solved using aiiy algorithm valid for the time-dependent Schrodinger 
equation (see e.g. Refs. j22|, 23, 24, 25, [2^, 13]). However, contrary to time-dependent 
models, it is obtained without semiclassical approximation: the projectile-target coordinate 
components b and Z are quantal variables in DEA. This advantage over time- dependent 
techniques allows taking into account interferences between solutions obtained at different bs. 
The DEA reproduces various breakup observables quite accurately for collisions of loosely- 



bound projectiles on both light and heavy targets [29|, [30 



The second step in the usual eikonal model is to assume the collision to occur during a 
very brief time and to consider the internal coordinates of the projectile to be frozen while 



the reaction takes place [17[. This second assumption, known as the adiabatic, or sudden, 
approximation leads to neglect the term Hq — Eq in Eq. fIlOp which then reads 



d - 

ihv—^(b, Z,r) 



VpT{R,r)^{b, Z,r) 



(11) 



In these notations, the asymptotic condition ([6]) becomes 



$(6,Z,r) — . e^''Hi^(«-^)l0/o.on.o(^o,r). 



The solution of Eq. (fTTi) exhibits the well-known eikonal expression 1^ 



(12) 



\1'(6, Z, r) = exp 



-— / VpT{b,Z',r)dZ' 



(13) 



This expression is only valid for short-range potentials. The Coulomb interaction requires a 
special treatment that is detailed in the next section. Let us point out that this treatment 
allows taking properly account of the projectile-target Rutherford scattering. The Coulomb 
distortion in Eq. ( |T2l) is therefore simulated in the phase of Eq. ( |T3l) . After the collision, the 
whole information about the change in the projectile wave function is thus contained in the 
phase shift x that reads 



X{b,s) = -^ [ VpT{R,r)dZ. 



(14) 



Due to translation invariance, this eikonal phase depends only on the transverse components 
b of the projectile-target coordinate R and s of the core-fragment coordinate r. 
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B. Coulomb correction to the eikonal model 



The eikonal model gives excellent results for nuclear-dominated reactions 17, 2^. How- 
ever, it suffers from two divergence problems when the Coulomb interaction becomes signif- 
icant. The first is the well-known logarithmic divergence of the eikonal phase describing the 
Coulomb elastic scattering 1^, 17, Tq]. The second is caused by the adiabatic approximation 



used in the eikonal treatment of the Coulomb breakup [IT]. To explain this, let us divide 
the eikonal phase f|T^ into its nuclear, and Coulomb contributions 

s) = s) + x^ib, s) + x'pAb). (15) 

The Coulomb term x'" for a one-neutron halo nucleus reads (the extension to the case of a 



charged fragment is immediate) j29|, |35 



mp 



= In 1-2^-- + ^- , 17 
\ mp rrip I 

b denotes a unit vector along the transverse coordinate b. In Eq. (1161) . we subtract the term 
1/R corresponding to a Coulomb interaction between the projectile center of mass and the 
target. The phase x'" therefore corresponds to the Coulomb tidal force that contributes to 
the breakup. Moreover, this subtraction leads to a faster decrease of the potential at large 
distances, which enables us to obtain the analytic expression f|T71) . This is compensated by 
the addition of the elastic Coulomb phase Xpt 

f • (18) 

This phase describes the Rutherford scattering between the projectile and the target. The 
integral is truncated, for it otherwise diverges (note that the integral in Eq. f[T7|) does 
not diverge and therefore does not require the same treatment). This truncation basically 
corresponds to Glauber's screened Coulomb potential [l6| . Other truncation techniques [l6| 
and other ways to deal with this divergence [l8| exist. All lead to the same expression of the 
elastic Coulomb phase but for an additional constant phase that does not affect the cross 
sections 1^. The truncation considered in Eq. (ITSl) leads to 

X^pAb) ^ 2r^\n-^. (19) 

This elastic Coulomb phase correctly reproduces Rutherford scatterin g, in dicating that the 



first of the two aforementioned divergences can be easily corrected |16l. Il7l. |18[|. The nuclear 
term x^ is then by definition the difference between the eikonal phase (fT^ and the Coulomb 
contributions (fT7|) and (fT9|) . 

In addition to the divergence in elastic scattering, the Coulomb interaction is responsible 
for a divergence in breakup. The aim of the present paper is to analyse a way to correct 
this divergence. It is due to the slow decrease of x*" in b. Indeed, when expanded in powers 
of x*"; the exponential of the Coulomb eikonal phase reads 

e'^" = l + ^X^-i(x^)' + ---, (20) 
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where the exphcit dependence on the coordinates has been omitted for clarity. When in- 
tegrated over b in the calculation of the cross sections (see Sec. Ill C\} . the 1/6 asymptotic 
behavior of the first-order term ix'" will lead to divergence. 

This divergence problem arises from the incompatibility between the infinite range of the 
Coulomb interaction and the adiabatic, or sudden, approximation: no short collision time 
can be assumed if the Coulomb interaction dominates. Renouncing the use of the adiabatic 
approximation solves this divergence: the DEA, which corresponds to the eikonal model 



without this approximation [see Eq. ffTOl) and Refs. |28|, l29j], does not diverge. The excellent 
results obtained within the DEA for collisions of loosely-bound projectiles on both light and 
heavy targets 2^, sO] confirm that, when dynamical effects are considered, both nuclear and 



Coulomb interactions can be properly taken into account on the same footing. 

To avoid this divergence, a cutoff at large b could be made. In Ref. [38i], Abu-Ibrahim and 
Suzuki proposed to limit the values of b to be considered in the cross-section calculations at 

hv , . 

This cutoff is obtained by requiring the characteristic time of internal excitation h/\Eo\ to 
be shorter than the collision time b/v. The factor of two is proposed as a qualitative guide. 
However this treatment is rather artificial and not very satisfactory [35'i]. 

Alternatively, it has been proposed by Margueron, Bonaccorso, and Brink js^], and de- 
veloped by Abu- Ibrahim and Suzuki [35|, to replace the first-order term ix'^ in Eq. (jgOl), 
which leads to the divergence, by the first-order term of the perturbation theory ix^^ |36| 



oo 



(22) 

where u = {E — EQ)/h, with E the c-f relative energy after dissociation. Since no adia- 
batic approximation is made in perturbation theory, this term does not diverge. When the 
adiabatic approximation is applied to Eq. (1221) . i.e. when u is set to 0, one recovers exactly 
the Coulomb eikonal phase f|T6l) . This suggests that without adiabatic approximation the 
first-order term in Eq. fl20l) would be ix^*^ fl22|) . Furthermore, a simple analytic expres- 
sion is available for each of the Coulomb multipoles in the far-field approximation, i.e. for 



m /r /m p < R [39j . The idea of the correction is therefore to replace the exponential of the 



eikonal phase according to 

e^x ^ e^x^ Lix"" _ ^^c _^ -^fo^, ^x%t _ (23) 



With this Coulomb correction, the breakup of halo-nuclei can be described within the 
eikonal model taking on (nearly) the same footing both Coulomb and nuclear interactions at 
all orders. This correction can also be seen as an inexpensive way to introduce higher-order 
effects and nuclear interactions in the first-order perturbation theory. 

In this work, we analyse the validity of this CCE model by comparing results obtained 
with the correction fl23p to results of the DEA. The latter is chosen as reference calculation, 
since it does not make use of the adiabatic approximation that leads to the divergence in 



the eikonal description of breakup. It is also in good agreement with experiments [29|, [30 
Calculations performed in the usual eikonal model, and at the first-order of the perturbation 
theory will also be presented to emphasize the effects of the correction. We focus on the case 
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of ^^Be breakup. In that case, only the dipole term of the Coulomb interaction is significant 
We thus restrict the correction to that multipole. The perturbative correction then 
reads ISal 



mp V 



(24) 



where K„ are modified Bessel functions 37 



Of course, in other cases, like in Coulomb 



breakup, the quadrupole term may no longer be negligible [l4,l3Q|, it should then be included 
in the correction. 



C. Breakup cross sections 



To evaluate breakup cross sections within the CCE we proceed as explained in Ref. [29 
replacing the DEA breakup amplitude by 



,c 



FO 



(25) 



where ai is the Coulomb phase shift ^]. The breakup amplitudes for the usual eikonal 
model are obtained in the same way but without the correction. 

In the following, we consider two breakup observables. The first is the breakup cross 
section as a function of the c-f relative energy E after dissociation [see Eq. (52) of Ref. j29| ] 



da 



bu 



dE t9k 2jo + 



1 fOO 



(26) 



mo Ijm 



This energy distribution is the observable usually measured in breakup experiments [Uj, \12 
It corresponds to an incoherent sum of breakup probabilities computed at each b 



h^k 2jo + 



(mo) 
kljm 



(27) 



mo Ijm 



The second breakup observable is the parallel-momentum distribution [see Eq. (53) of 



Ref. 29 



dcT] 



bu 



dk\\ 



2j 



OTT X ^ 

mo 



bdb 



dk 



^(//m - vv\3m)Yr-^{e,, o)s!;;°2{b) 



(28) 



where dk = aiccos{k / k\\) is the colatitude of the c-f relative wavevector k after breakup. 
Contrary to the energy distribution, the parallel-momentum distribution corresponds to a 
coherent sum of breakup amplitudes. This observable is therefore sensitive to interferences 
between different partial waves. Consequently, it constitutes a particularly severe test for 
reaction models. 
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III. NUMERICAL ASPECTS 



For these calculations, we use the same description of ^^Be as in Ref. 4l|. The halo 
nucleus is seen as a neutron loosely bound to a ^°Be core in its 0^ ground state. The ^°Be-n 
interaction is simulated by a Woods-Saxon potential plus a spin-orbit coupling term (see 
Sec. IV A of Ref. 4l|)- The potential is adjusted to reproduce the first three levels of the 
^^Be spectrum. The ^ ground state is seen as a lsl/2 state, while the | excited state 
is described by a Opl/2 state. This well-known shell inversion is obtained by considering a 
parity-dependent depth of the central term of the potential. The | resonance at 1.274 MeV 
above the one-neutron separation threshold is simulated in the d'b/2 partial wave. 

The interaction between the projectile components and the target are simulated by opti- 
cal potentials chosen in the literature. In our calculations, we use the same potentials as in 



Refs. |27l . |41| . As suggested in Ref. [42|, the Be-Pb potential is scaled from a parametri- 



sation of Bonin et al. |43| that describes elastic scattering of 699 MeV a particles on lead 
[potential (1) in Table III of Ref. j27i|]. For the ^°Be-C interaction, we use the potential 
developed by Al-Khalili, Tostevin, and Brooke, which reproduces the elastic scattering of 
^°Be on C at 59.4 MeV/nucleon Q (potential ATB in Table III of Ref. I4l|). In both cases, 
we neglect the possible energy dependence of the potential. We model the neutron-target 



interaction with the Becchetti and Greenlees parametrisation [45 



To evaluate the breakup amplitude (!25|) within the CCE or the usual eikonal model, we 
need to compute the eikonal phase (fT5l) . The nuclear part is evaluated numerically, while 
the Coulomb part is obtained from its analytic expression (fT7|) . The numerical integral over 
Z is performed on a uniform mesh from ^min = ~20 fm up to 2^max = 20 fm with step 
AZ = 1 fm. The corrected phase (1231) is then numerically expanded into multipoles of rank 
A. We use a Gauss quadrature on the unit sphere similar to the one considered to solve the 
time-dependent Schrodinger equation in Ref. 23|. The number of points along the colatitude 
is set to Ng = 12, and the number of points along the azimuthal angle is = 30. Unless 
otherwise stated, we perform all calculations with multipoles up to Amax = 12. 

The eigenfunctions of the projectile Hamiltonian Hq ([T]) are computed numerically with 
the Numerov method using 1000 radial points equally spaced from r = up to r = 100 fm. 
The same grid is used to compute the radial integral in Eq. fl25|) . For Coulomb (nuclear) 
breakup, the integrals over b appearing in Eqs. fl2^ and fl25]) are performed numerically 
from = up to 6 = 300 (100) fm with a step Ab = 0.5 (0.25) fm. 

The DEA Schrodinger equation flTUl) is solved using the numerical technique detailed in 



Ref. [27|. In this technique, the projectile internal wave function is expanded upon a three- 
dimensional spherical mesh. The size of the mesh required for the calculation varies with 
the projectile-target interaction. For Coulomb (nuclear)-dominated reactions, the angular 
grid contains up to Ng = 8 (12) points along the colatitude 9, and N^, = 15 (23) points 
along the azimuthal angle (f. This corresponds to an angular basis that includes all possible 
spherical harmonics up to / = 7 (11). The radial variable r is discretized on a quasiuniform 
mesh that contains A^^ = 800 (600) points and extends up to r^r^ = 800 (600) fm. The time 
propagation is performed with a second-order approximation of the evolution operator. It is 
started at tin = —20 (10) h/MeV with the projectile in its initial bound state, and is stopped 
at tout = 20 (10) h/MeV (t = corresponds to the time of closest approach). The time step 
is set to At = 0.02 h/MeV in both Coulomb and nuclear cases. 

The evolution calculations are performed for different values of b. These values range from 
up to 300 (100) fm with a step Ab varying from 0.5 (0.25) fm to 5.0 (2.0) fm, depending 
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b (fill) 

FIG. 2: Breakup probabilities as a function of transverse coordinate b for ^^Be impinging on 208p|-, 
at 69 MeV/nucleon. Three energies E are shown: 0.5 MeV, 1.274 MeV, and 3.0 MeV. The results 
are obtained within DEA (full lines), CCE (dotted lines), usual eikonal approximation (dashed 
lines), and first-order perturbation theory (dash-dotted lines). The upper part displays the values 
at small b, while the lower part focuses on the asymptotic region. 



on b. The integrals over b are performed numerically. 



IV. BREAKUP OF "BE ON PB AT 69 MEV/NUCLEON 

We first consider the breakup of ^^Be on lead at 69 MeV/nucleon, which corresponds to 



the experiment of Fukuda et al. at RIKEN [12| . These data are fairly well reproduced by the 
DEA H, that we use as reference calculation. Since we focus on the comparison of models, 
we do not display Fukuda's measurements. A comparison with experiment would indeed 
require a convolution of our results, which would hinder the comparison between theories. 

In Fig. [21 we compare the breakup probability fl271) obtained with the DEA (full lines), 
the CCE (dotted lines), the usual eikonal model (Eik., dashed lines), and the first-order 
perturbation theory (FO, dash-dotted line). They are depicted as a function of the transverse 
coordinate b for three ^°Be-n relative energies: E = 0.5 MeV, 1.274 MeV (i.e. the | 
resonance energy), and 3.0 MeV. The upper part of Fig. [2] displays the values at small b, 
while the lower part, in a semilogarithmic scale, focuses on the asymptotic region. 
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Over the whole range in b, the CCE results are close to the DEA ones, and this at all 
energies. This good agreement suggests the Coulomb correction to be valid for simulating the 
breakup of loosely-bound nuclei on heavy targets. In particular, the CCE is superimposed 
to the DEA results in the asymptotic region. Obviously, the first-order perturbation theory 
efficiently corrects the erroneous 1/6 asymptotic behavior of the usual eikonal model. 

At small b, the agreement between the CCE and DEA seems slightly less good. In par- 
ticular, at small energy, the corrected eikonal model overestimates the reference calculation. 
This is due to the far-field approximation used in the first-order perturbation correction. This 
approximation provides a convenient analytical expression (1241) of the phase x^'^- However, 
it is incorrect at small b: it diverges at 6 = 0. Nevertheless, in spite of that divergence, 
the CCE remains close to the DEA. This illustrates that the CCE can also be seen as a 
way to include nuclear interactions within the first-order perturbation theory, and correct 
its ill-behavior at small b. 

The breakup cross section fl2^ computed with the four approximations is displayed in 
Fig. El^a) as a function of the ^''Be-n relative energy E after dissociation. Contributions of 
the s, p, and d partial waves are shown separately in Fig. El^b). The small bump at about 
1.25 MeV is due to the resonance in the d5/2 partial wave. The CCE cross section (dotted 
line) is nearly superimposed on the DEA one (full line). Only at low energy is the CCE 
slightly larger than the reference calculation. As mentioned earlier this effect is due to the 
use of the far- field approximation to derive the perturbative correction x^"^- 

Interestingly, the agreement between CCE and DEA is better for the total cross section 
than for each partial-wave contribution: The CCE p contribution is larger than the DEA 
one, while the CCE s and d contributions are smaller than the DEA ones. We interpret this 
as a lack of couplings in the continuum in the CCE. In the DEA, these couplings depopulate 
the p waves towards the s and d ones without modifying the total cross section [40]. The 
differences between CCE and DEA partial-wave contributions suggest that this mechanism 
is hindered in the former. 

The wrong asymptotic behavior of the Coulomb eikonal phase f|T7|) leads to a divergence 
in the calculation of the breakup cross sections. To evaluate the energy distribution within 
the usual eikonal model one needs to resort to a cutoff at large b. The cutoff proposed in 
Ref. [see also Eq. fl?Il) ] gives here 6max = 71 fm. The corresponding cross section is 
displayed in Fig. [3](a) with a dashed line. Its energy dependence is strongly different from 
that of the reference calculation: it is too small at low energy and too large at high energy. 
The p contribution, which includes the diverging term of the Coulomb eikonal phase (fT7|) . is 
responsible for that ill-behavior. Contrarily, the s and d contributions are superimposed on 
those of the CCE. The use of the Coulomb correction therefore significantly improves the 
eikonal model when considering collisions with heavy targets. 

The cross section obtained within the first-order perturbation theory is shown in dot- 
dashed line. The nuclear interactions between the projectile and the target are described by 
a mere cutoff at 6min = 15 fm. This value has been chosen to fit the DEA energy distribution 
in the region of the maximum. Here again, the shape of the cross section is very different 
from that of the reference calculation. However, contrary to the usual eikonal model, it 
decreases too quickly with the energy. Moreover, since only the dipole term of the Coulomb 
interaction is considered, only the p wave is reached from the s ground state, whereas s and 
d waves are significantly populated through nuclear interactions and higher-order effects. 
Note that a smaller cutoff 6min, in better agreement with the usual choice that corresponds 
to the sum of the projectile and target radii, does not improve the agreement. 
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E (MeV) 

FIG. 3: (a) Breakup cross sections for ^^Be impinging on ^'^^Pb at 69 MeV/nucleon as a function 
of the ^"Be-n relative energy E. The results are obtained within the DEA, the CCE, the usual 
eikonal approximation with upper cutoff 6jnax = 71 fm, and the first-order perturbation theory 
with lower cutoff 6min = 15 fm. (b) Contributions of the s, p, and d partial waves. 

We now consider the parallel- momentum distribution [see Eq. (128!) ]. This breakup observ- 
able is more sensitive to interferences and therefore constitutes a more severe test than the 
energy distribution. The parallel-momentum distribution computed within the four models 
is displayed in Fig. HI 

As in the previous cases, the CCE is in excellent agreement with the DEA in both 
magnitude and shape. We simply note that the former is slightly less asymmetric than the 
latter, which is probably a signature of the lack of couplings in the continuum mentioned 
earlier. On the contrary, both the usual eikonal model and the first-order perturbation 
theory lead to rather poor estimates of the momentum distribution. First, they lead to an 
erroneous magnitude of the cross section. The usual eikonal model gives too large a parallel- 
momentum distribution. This is related to the too slow decrease obtained for the energy 
distribution. On the contrary, the first-order perturbation gives too low a cross section; 
a defect due to the quick decrease in the energy distribution. Lowering the cutoff bmm to 
cure this problem would then lead to too large an energy distribution in the peak region. 
Second, none of these models exhibits the asymmetry observed in the DEA. This absence of 
asymmetry in parallel-momentum distributions of the breakup of loosely-bound projectiles is 
a well-known problem of the eikonal model [46j] . It is fortunate that the Coulomb correction, 
combining two approximations that lead to perfectly symmetric momentum distributions, 
restores the asymmetry observed experimentally and in dynamical calculations. 
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FIG. 4: Breakup cross sections for ^^Be impinging on ^''^Pb at 69 MeV/nucleon as a function of the 
^^Be-n relative parallel momentum Aiy. The figure displays the results obtained within the DEA, 
the CCE, the usual eikonal approximation with an upper cutoff 6max = 71 fm, and the first-order 
perturbation theory with a lower cutoff ftmin = 15 fm. 
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FIG. 5: Convergence of the multipole expansion in Amax of the CCE illustrated on the parallel- 
momentum distribution computed for ^^Be impinging on 208p^ gg MeV/nucleon. 

Fig. illustrates the convergence of the CCE with regard to the number of multipoles 
considered in the breakup computation. The parallel-momentum distributions obtained with 
maximum multipolarities Amax = 4, 8, and 12 are displayed. Although all three calculations 
are close to one another, Amax = 4 has not yet converged: there remains some 4% difference 
with the other two at the maximum. On the contrary, the difference between Amax = 8 
and 12 is insignificant (about 0.5%). This shows the necessity to include a large number of 
partial waves in dynamical calculations. Note that other breakup observables converge with 
a lower number of multipoles. In particular, the energy distribution requires only Amax = 4 
to reach satisfactory convergence. 

These results confirm the ability of the Coulomb correction to reliably reproduce breakup 
observables for collisions of loosely-bound projectiles on heavy targets. It reproduces dy- 
namical calculations with an accuracy that is unreachable within the usual eikonal model 
or the first-order perturbation theory, on which it is based. 
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V. BREAKUP OF ^^BE ON C AT 67 MEV/NUCLEON 

To complete this analysis of the Coulomb correction, we investigate its effect in nuclear 
induced breakup. The usual eikonal description of such reactions is known to give excellent 
results (itI . 29| . The Coulomb interaction between the projectile and the target plays then 



a minor role and we expect the correction fl2^ to have much less influence than in the 
Coulomb breakup case. 

For this analysis, we consider the breakup of ^^Be on a carbon target at 67 MeV/nucleon, 
which corresponds to the experiment of Fukuda et al. T3]. The DEA is in excellent agree- 



ment with Fukuda's data [29j, and therefore constitutes our reference calculation. For the 
same reasons as in the previous section, we do not compare directly our calculations with 
experiment. 

Fig. [6] displays the breakup probability (1271) obtained at three energies E = 0.5, 1.274, 
and 3.0 MeV within the DEA (full lines), the CCE (dotted lines), and the usual eikonal 
model (dashed lines). Since this reaction is nuclear dominated, we no longer display the 
result of the first-order perturbation theory. The upper part of Fig. E] displays the breakup 
probability at small b, while the lower part emphasizes the asymptotic behavior of Pbu in a 
semilogarithmic plot. 

In this case, all three reaction models lead to similar results. This confirms the validity 
of the adiabatic approximation in the eikonal description of nuclear-dominated reactions. 
The difference between the DEA and the other two models is indeed rather small. Only 
at E = 1.274 MeV, the energy of the I"*" resonance, does it become significant (up to 10% 



difference in the vicinity of the peak at 6 ~ 6 fm). This larger difference suggests stronger 
dynamical effects at the resonance. This is not very surprising since the presence of that 
resonance strongly increases the breakup process [4lj |. 

Up to 6 = 20 fm, the usual eikonal model and the CCE remain very close to one another, 
confirming the small role played by the Coulomb interaction in the dissociation. At larger 
b, where only Coulomb is significant, we observe the 1/6 behavior of the usual eikonal 
model. This ill-behavior is corrected using the CCE, whose breakup probabilities are nearly 
superimposed on the DEA ones in the asymptotic region. However, since this correction 
affects breakup probabilities at two or three orders of magnitude below the maximum, we 
do not expect it to significantly influence breakup observables. 

The breakup cross sections computed within the three models are plotted as functions 
of the energy E in Fig. [71 The contributions to the total cross section of the partial waves 
s, p, and d are shown as well. The large peak at about 1.25 MeV is the signature of the 
significant enhancement of the breakup process by the d5/2 resonance. As suggested by the 
previous result, all three models lead to very similar cross sections. This similarity is also 
observed in the partial-wave contributions. The couplings in the continuum that depopulate 
one partial wave toward others, as observed in Coulomb breakup (see Fig. [3] and Ref. (4o|). 
are thus much smaller in nuclear-induced breakup. 

As in Fig. [HI the difference between the DEA and the other two models is rather small. 
The DEA is about 6% in average larger than the eikonal model. Note that this difference 
reaches 8% at the resonance energy, which is consistent with the difference observed in 
Fig. [6](a). The usual eikonal and the CCE lie even closer to one another. The relative 
difference between them in the total cross section does not exceed 3%. Even in the p partial 
wave, where the Coulomb correction is performed, no significant difference is observed. This 
confirms that the correction of the eikonal model is not necessary for nuclear-dominated 
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FIG. 6: Breakup probabilities as a function of transverse coordinate b for ^^Be impinging on ^^C at 
67 MeV/nucleon. Three energies E are shown: 0.5 MeV, 1.274 MeV, and 3.0 MeV. The results are 
obtained within the DEA (full lines), CCE (dotted lines), and usual eikonal (dashed lines) models. 
The upper part displays the values at small b, while the lower part emphasizes the behavior in the 
asymptotic region. 



reactions due to the small role played by the Coulomb interaction. The cutoff in b proposed 
in Ref. [sl] is therefore sufficient. 

The parallel-momentum distributions obtained with the three models are displayed in 
Fig. [HI As already mentioned, this observable is a more severe test for reaction models than 
the energy distribution. We observe significant differences between the DEA and the other 
two models. As in the case of Coulomb breakup, the DEA leads to an asymmetric parallel- 
momentum distribution: The DEA distribution is shifted toward neg ative fcy and presents 



a more developed tail on the negative fcy side, as observed in Ref. [46 . 

As for the previous observable, the CCE and usual eikonal models lead to very similar 
parallel-momentum distributions. These distributions are symmetric. As mentioned earlier, 
this symmetry is due to the lack of dynamical effects in the eikonal description of reactions. 
Contrary to the Coulomb case, the correction fl25]) is not able to restore this asymmetry. 
It indicates that these dynamical effects result from the nuclear interactions between the 
projectile and the target. 

The convergence of the CCE model with the number of multipoles is illustrated in Fig. [9] 
for the parallel-momentum distribution. The CCE distributions computed with Amax = 4- 
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FIG. 7: Breakup cross sections for ^^Be impinging on ^^C at 67 MeV/nucleon as a function of the 
^°Be-n relative energy E. Results are obtained within the DEA, the CCE, and the usual eikonal 
model with an upper cutoff 6max = 70 fm. Contributions of the s, p, and d partial waves are shown 
as well. 
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FIG. 8: Breakup cross sections for ^^Be impinging on ^^C at 67 MeV/nucleon as a function of the 
^'^Be-n relative parallel momentum /cy. Results are obtained within the DEA, the CCE, and the 
usual eikonal approximation with an upper cutoff 6max = 70 fm. 



12 are displayed. The convergence is much slower than for Coulomb-dominated breakup 
(see Fig. [5]). The relative difference between Amax = 10 and Amax = 12 is indeed about 
3% at the maximum. This is due to the rapid variation of the nuclear potential with the 
projectile-target coordinates. It confirms the need of a larger number of partial waves in 
the dynamical calculation of nuclear-dominated dissociation. Note that the convergence is 
faster for the energy distribution. For that observable, an acceptable convergence is reached 

-^max ~ 6. 
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FIG. 9: Convergence in Amax of the CCE illustrated on the parallel- momentum distribution for the 
breakup of ^^Be on ^^C at 67 MeV/nucleon. 



VI. CONCLUSION 

The eikonal description of reactions is a useful tool to simulate breakup and stripping 



reactions on light targets at intermediate and high energies |16l . Il7l . |29| . This model is inter- 
esting because of its relative simplicity in implementation and interpretation with respect 
to other elaborate models, like CDCC or DEA. Unfortunately, it suffers from a divergence 
problem associated with the treatment of the Coulomb interaction between the projectile 
and the target. This divergence is due to the incompatibility of the adiabatic, or sudden, ap- 
proximation which is made in the usual eikonal model, and the infinite range of the Coulomb 
interaction. One way to cure this problem is not to make this adiabatic approximation. This 
leads to the DEA 28, 2^. However, like other elaborate models, the DEA is computationally 
expensive. Another way to solve this problem is to substitute the diverging Coulomb phase 
at the first-order of the eikonal model by the corresponding first-order of the perturbation 



theory [3J, |35 



In this work, we study the validity of this Coulomb correction by comparing it to the 
DEA, which does not present the divergence problem of the usual eikonal model. The 
chosen test cases are the dissociation of ^^Be on Pb and C at about 70 MeV/nucleon. These 



correspond to RIKEN experiments [ll|, Il2[ that are well reproduced by the DEA [2t 

In the case of the Coulomb breakup, the CCE gives results in excellent agreement with the 
DEA. The combination of the eikonal model with the first-order perturbation theory indeed 
solves the divergence problem due to the Coulomb interaction. Moreover, it correctly takes 
into account the nuclear interaction between the projectile and target. The breakup observ- 
ables (energy and parallel-momentum distributions) obtained within the DEA are accurately 
reproduced using the CCE. This agreement is obtained while both CCE ingredients — usual 
eikonal and first-order perturbation — fail to describe the reaction. First they both require 
a rather arbitrary upper or lower cutoff in b in order not to diverge. Second they do not re- 
produce the shape of the breakup cross sections. In particular the CCE gives an asymmetric 
parallel-momentum distribution, in agreement with the dynamical calculation. Contrarily, 
both the usual eikonal and the perturbative models lead to perfectly symmetric distributions. 
This suggests that CCE restores dynamical effects that are missing in its ingredients. 

The Coulomb correction has much less effect on the nuclear-dominated breakup. This 
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was expected because of the much smaller influence of the Coulomb interaction in reactions 
involving light targets. This result indicates that in this case the correction is not essential. 
It also implies that the CCE suffers the same lack of dynamical effects as the usual eikonal 
model in nuclear dominated reactions. 

The CCE successfully combines the positive aspects of both the eikonal model and the 
first-order perturbation theory. It allows describing accurately the nuclear interaction while 
correctly reproducing Coulomb-induced effects. Moreover the CCE restores some of the 
dynamical effects, which are totally absent in other simple models. It therefore provides 
a reliable description of the breakup of loosely-bound projectiles at intermediate and high 
energies. Its simplicity in use and interpretation suggests it as a competitive alternative to 
more elaborate models to describe the breakup of Borromean nuclei. 
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